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Abstract 

I present an asymptotic formula for the Takeuchi numbers T„. In particular, I 
give compelling numerical evidence and present a heuristic argument showing that 

Tn^CTBne^v\w{nf 

as n tends to infinity, where are the Bell numbers, W{n) is Lambert's W function, 
and Ct = 2.239 ... is a constant. Moreover, I show that the method presented here 
can be generalized to derive conjectures for related problems. 
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1 Introduction 



In a paper entitled "Textbook Examples of Recursion," Donald E. Knutli discusses re- 
currence equations related to the properties of recursive programs [[Knutli 1991|| , among 
them Takeuchi's function [ [I'akeuchi 1978| , [I'akeuchi 1979 | 



t{x, y,z) = if X < y then y else t(t{x — 1, y, z),t{y — 1, z, x),t{z — 1, x, y)) . (1) 

Let T{x,y,z) denote the number of times the else clause is invoked when t{x,y,z) is 
evaluated recursively.^] For non-negative integers n, the Takeuchi numbers T„ are defined 
as Tn = T{n, 0,n + 1). The first few values of T„ for n = 0, 1, 2, . . . are 

0, 1, 4, 14, 53, 223, 1034, 5221, 28437, 165859, .... (2) 

Knuth gives the recurrence 

r /n + A;\ fn + k\]^ f2k\ 1 



fc=0 



and deduces a functional equation for the generating function T{z) = Z^^o^"-^"- 

m ^ ^ . ^<^n.ci.y> . (4) 

I — z V 1 — 4z 



where 



2n\ 1 



c'(^) = t:(i-v'i^) = E L h-^ (5) 



22 ^o\n Jn + 1 



is the generating function for the Catalan numbers C„ = Lastly, he gives asymp- 

totically vahd bounds for T„, 

gnlogn-„loglogn-n < 7; < gnlogn-n+logn ^jj sufficiently large U, (6) 

and poses obtaining further information about the asymptotic properties of T„ as an open 
problem. 



^Note that it is the recursive evaluation of t{x, y, z) rather than the actual value of t{x, y, z) that is of 
interest. See Knuth's paper for an explicit expression of t(a;, y, z). 



In this paper I give arguments leading to two conjectures about the asymptotic be- 
haviour of the Takeuchi numbers T„. In Section 2, 1 present an explicit asymptotic formula 
for T„ which improves upon the bounds (P) based on numerical evidence and a heuristic 
argument. For this, I briefly discuss the related asymptotic behavior of the Bell numbers 
and give an argument based on a numerical observation which leads directly to an explicit 
asymptotic formula for T„ as n tends to infinity. The formula, as described in Conjecture 
1, is exact up to 0((logn/n)^) and contains a constant Ct which is numerically deter- 
mined to 25 significant digits. Section 3 presents a heuristic analytic argument which 
gives the asymptotic behavior up to o(l) and enables one to identify the constant Ct in 
terms of an explicit expression, as stated in Conjecture 2. In the final section I show 
that the method developed in Section 3 can give insight into the asymptotic behavior of 
a larger class of problems. 

I conclude this introduction by briefly discussing the structure of the recurrence (|^) 
and the related functional equation (H). It is clear from the asymptotic bounds (H) for T„, 
that the generating function T(z), defined as a formal power series, does not converge, and 
is therefore at best only an asymptotic expansion to an actual solution of the functional 
equation. This is also evident from the structure of the functional equation. This structure 
becomes clearer upon a change of variables, 

nz) = -T{z - z') - ^ , (7) 

z [1 — z)[l — z + z^) 

where one sees directly that the transformation involved is g{z) = z — z^, which is only 
marginally contracting at its fixed point z = 0. While functional equations with a trans- 
formation g{z) that has an expanding or contracting fixed point (|5''(0)| ^ 1) are very well 
understood ||Kuczma 1990|| , it is precisely the fact that |5''(0)| = 1 which is at the root of 
the underlying difficulty of the problem discussed in this paper. 
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2 Numerical Observations 



Our starting point is Knuth's observation [[Knuth 1991|| that for n > the Bell numbers 



Bri are a lower bound to T„. Here, Br, is defined as 



k=0 



The asymptotics of -B„ is discussed in great detail by de Bruijn i DeBruijn 1961 ]. In 



[Moser 1955|| one finds a systematic way of generating higher order terms in the asymptotic 



expansion by means of a contour integral representation using the well-known fact that 

Y: B^- = exp(e^ - 1) . (9) 



n=o 



Alternatively, one can also expand the right-hand side of in z to get 

Bn=-T.^^ (10) 

which can then be evaluated asymptotically by use of the Euler-MacLaurin formula. Ei- 
ther way, in the course of the computation of the asymptotics of i?„ it turns out that a 
convenient asymptotic scale is given in terms of W{n) rather than in terms of n, where 
W{x) is Lambert's W function, which is defined as the real solution of 

W{x)expW{x) = X . (11) 

The sum ([T0|) is dominated by terms around m = e^^"\ and one can easily calculate that, 
written in terms oi w = W{n), the Bell numbers B^ behave asymptotically as 

logB„ = e-lu.^ - u. + 1) - 1 log(l + u.) - 1 - ^^'^l^^ll^ e- (12) 

w{2w^ + 12w^ + 29w^ + A0W + 36) -3w^ 



48(l + ^^;)6 

and it is straightforward to calculate additional terms. It is more conventional to state 
this formula with the exponentials e^^ replaced by (n/w)'', but this obscures the fact that 
the asymptotic expansion is obtained in terms of w rather than n. 



Having rather explicit control over this lower bound, it is natural to now try to compare 
Bn and T„ more closely. There is a principal difficulty coming from the fact that the 
asymptotic scale presumably also involves w — W{n), which grows more slowly than 
logn. (As an example, Vr(lOOO) « 5.2496 and 1^(10000) a; 7.2318.) Thus, one would 
expect that a direct numerical investigation of T„/_B„ is not very insightful, due to the 
presence of slowly varying correction terms of unknown form. 
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Figure 1: T„+i/T„ — B^/ Bn-i plotted versus 1/n for n < 1000 
However, as Figure 1 shows, if one compares the growth rates T„/T„_i and B^/ B^-i 
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Figure 2: Tn+i/[Bnexp{^w'^ + w)] plotted versus 1/n for n < 1000. The horizontal line 
is at Ct = 2.2394331040 . . .. 

instead, one is led to observe the surprisingly simple relationship 



lim 



T 



Bn-l 



1 . 



(13) 



In fact, the left hand side approaches 1 rather quickly, 

-1 + 0(6-^") 



(14) 



This (unproven) numerical observation leads to a straightforward derivation of an asymp- 
totic formula. From (|12D it follows easily that = e""" + 0(6"^"'). Now one takes 
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logarithms and sums up successively, from whence it follows that 

logT„+i = log5„ + ]^w^ + w + 0(1) . (15) 

Now that I have guessed the leading asymptotic form, 1 can again resort to numerical 
work to try to improve upon it. In fact, numerically it appears that the convergence is 
even better than expected due to a chance cancellation of higher order correction terms. 
Figure 2 shows that 

Tn+i = Ct Bn expi^w^ + W + ©(e'^-)) , (16) 

and from the first 1000 series terms I am able to deduce by iterative application of standard 
series extrapolation methods that 

Ct = 2.23943 31040 05260 73175 4785 (1) . (17) 

Using the known asymptotic form of the Bell numbers, I can now give an explicit 
asymptotic expression for T„ in terms of w = W{n) alone, as stated in the following 
conjecture. 

Conjecture 1 As n tends to infinity, one has 

\ogTn = e'^iw^ -w + l) + ^w^ -hog{l + w)+ (18) 



24(1 + wy 

Here w = W{n) and the constant Ct is some positive real number. 



Numerically, logC^ - 1 = -0.19377 72447 31916 75890 1157 (1). Of course it would 
be desirable to find an analytic expression for this number. In the next section I shall 
present a heuristic argument giving such an expression. 

Dropping one correction term and comparing with the asymptotic expression (|12D for 
the Bell numbers. Conjecture 1 implies the nice formula given in the abstract, 

T„~CT5„expi^y(r^)^ (19) 



3 Analytic Results 

In view of the previous section it seems promising to exploit the apparent affinity between 
Takeuchi numbers T„ and Bell numbers Bn- Given a recurrence of the general form 

n 

an = Yl Cn,ka'n-k + K , (20) 
k=l 

I choose to write 

1 °° 1 1 °° 1 

an = -Y\ —;fm,n and bn = -Y] —rbn ■ (21) 

Inserting (^) into the recurrence (|20|) and shifting the summation index by one, I next 
equate terms to get 

n 

fm,n ^ ^ (^n,kfm—l,n—k ~l~ ) (^2) 

A:=l 

This Ansatz might seem less arbitrary when considering that in the case of Bell numbers 
it reduces to fm,n = In general, one observes that fm,n must be a polynomial in m of 
at most degree n, which I write as 



n 



fm,n = m''Y. (^n,km ^ . (23) 

fc=0 

If one further requires the coefficients Cn^k in the recurrence to be polynomials of degree 
k in n, it follows that dn,k are polynomials of degree k in n, which I write as 

k 

dn,k = n''Y.rk,in-^ . (24) 

1=0 



Combining equations (p3|) and (|23) gives 



n k n n—l 

= Y.Y. rk,in''~^m'^ = Y^ni-^Y. ri+kA^/^f ■ (25) 

k=0 1=0 1=0 k=0 

In order to get an idea about the asymptotic behavior of this double sum, I now replace the 
quotient n/m by a new variable v and consider the formal limit of taking the summation 
bounds to infinity, leading to 

oo oo 

Smiv) = ^m^Wi{v) with r,(f) = J2'^i+k,iv'' ■ (26) 

1=0 k=0 



Applying this method to Takeuchi numbers, one now inserts Cn,k — { ("n^i^^) ~ ("^n~^) } 
and bn — Y,k=i (T) fe+T- With this choice, ro{v) is trivially zero as Cn,k are polynomials in 
n of degree /c — 1, and one gets a rather interesting result for Z > 1. In fact, 

ri{v) = (27) 
r^{v) = e^"' (^2e'''- i(T;3 + ^;2 + 4^; + 2)e^) (28) 

r3(^;) =. e'^''' (—e'^'' -{v^ + 3v^ + 7v + 6)e^''+ (29) 
V 8 

+ -^(3v^ + 6v^ + 47v^ + 52v^ + UAv'^ + 7Av + 51)6"^ 
r^{v) = e'^^" (^-^e^^ + ^{v' + 5v^ + 12v + 12)e^^+ (30) 
+ -^(3v^ + 18v'' + 89t;'^ + 226t;-^ + 411v'^ + 406v + 195)e^" - 
^ (9^^^ + 27v^ + 315v^ + 603v^ + 3024^;^ + 



432 

+3384t;^ + 8757^;^ + 4707^;^ + 5484v + 772) 
Prom this, I conjecture 

ri{v) = {piMe'^'+piMe^'-'^'' + • • . + pi,i-,{v)e'') , (31) 

where Pi,k{v) are polynomials in v of degree 3k. (This pattern has been verified for I < 8.) 
For V large, this conjecture implies 

ri{v) ~ Xie-y+'- , (32) 

and with a little effort one can compute the next values of A; = pi^Q. One gets 

1 ^_347 28201 
Ao-0, Ai-1, A2-2, A3---, A4--— , A5-^^, (33) 

_ _ 3172987 _ 822813607 _ 2183235065857 

^ ~ ~ 216000 ' ^~ 93312000 ' ^~ 16003008000 ' ^ ' 

from whence I am led to conjecture that one can write A^ = A*n/[(?T' — 1)']^ where /in is 

integer. (Unfortunately, I did not find a way to compute the coefficients A/ in a closed 

form!) One would expect from this behavior that 

Sm{v) ~ e"^^"h{eym) , (35) 



where h{x) is given in some sense by 

oo 

h{x) = J2\ix'. (36) 

1=0 

I caution here that the series may be divergent and just vahd as an asymptotic expansion. 

Keeping in mind that the evidence for the existence of h{x) is rather sketchy, I never- 
theless proceed under the assumption that ioi n » m » 1 one can write 

fm,n ~ m"e^("/"^)'/i(e"/"^/m) . (37) 

This now enables a heuristic computation of the Takeuchi numbers T„. I approximate 

1 °° 1 1 1 

Tn — - 22 — \fm,n ~ - /J — r/m,n (38) 



e , . ml 

mKimmax (n) 



^e^("/")'/i(e"/"7m) , (39) 



where in the last step it is assumed that n » m^a-r(n) >> 1. This sum is indeed 

dominated around ^moa;(^) ~ e"", where the argument of h simplifies to 1. A careful 
asymptotic analysis of 



^" = 7 E ^e^("/"^)'/i(e"/"^/m) (40) 

gives 



e / . ml 



logr„ = e"'(^/;2 - + 1) + - 1 log(l + w) + /io - 1 (41) 

^(12^5 + 2Aw^ + 36w=^ + 58^2 + 29w - 10) _^ 

^ 2A{w + 1)3 ^ 

^ (w+l)(/^f + /^2) + (2w^ + W + 2)/.i ^_^ ^ ^ 

where one has expanded h{x) around x = 1 as log h{x) — -\- hi{x — 1) -\- h2{x — 1)^/2 + 

o{{x-in 

As long as the corrections made on passing from T„ to Tn are small enough, it follows 
easily from this that asymptotically 

r„ ~ By^^'h{l) , (42) 
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and one can identify the constant Ct from equation (|T^) with h{l). Provided the series 
expansion of h{x) = J^kLo^k^'' converges at x = 1, I can thus conjecture an exphcit 
expression for the constant Ct, which in principle is computable. 

Conjecture 2 The constant Ct in Conjecture 1 is given by 

oo 

CT = hil) = Y,\k. (43) 

fc=0 

While the approximation of T„ by T„ may be correct up to 0{e^^), no choice of 
h{x) can match the next term in (^T]) with the expansion of T„. Thus, one also gets an 
indication of the size of the error made. 

It seems that a careful asymptotic evaluation of the fm,n promises to be a suitable 
way of providing rigorous proof for the asymptotics of the Takeuchi numbers. Of course 



one could also try to find a direct proof of our numerically observed equation (O) 



4 A Generalization 

In the derivation of the functional equation (|^) for the Takeuchi numbers T„, it is crucial 
that 

t["\'^^)^' = C(^f/VT^z, (44) 

as this identity allows the explicit summation of the terms in the recurrence (^. The 
identity used is a special case of the following nice identity 

This identity can be proved by inserting z = y/{l +?/)^"^^, which after expanding leads to 



k=0 

and 



^/(A + 1)A;\ z'' , , 



fc=0 



11 



I use this now as a motivation for the study of the family of recursions (with parameter 
A) 

Due to equation (|45|) one is able to derive a functional equation for the corresponding 
generating function A{z) = Z^J^Lo^"^"- 

A{z) = 1 + z^^A{z{l + y)) , z = y/{l + y)'+' . (49) 

For A = one recovers the recursion for the Bell numbers, and for A = 1 one has something 
which is at least "morally" related to the Takeuchi numbers. 

Inserting the Ansatz ( pi]) into (|i8|), one can easily repeat the analysis of the previous 
section. The result is now 

v4„~E„expA|iiy(n)' + l^(n) + c/(A)| (50) 

for any fixed value of A. Again, one have an identification of the kind d{X) = hx{l), where 
the first terms in the series expansion of hx{x) are 

hx{x) = i(A-l)s-^(2A2 + 18A-5)x2-^(33A^ + 90A2- 329A + 54)a;3 

^ -(52A^ - 520A^ + 4240A - 502)x^ + O(x^) , (51) 



960 

and one sees that the A;th coefficient is a polynomial in A of degree k (this has been 
verified up to A; = 7). I caution again that convergence of this series expansion is an open 
question. 

Finally, one can establish numerically the next term in the asymptotic expansion of 
An. For any fixed value of A, one finds 

log A„ = log i?„ + A + W7 + d{X) - ^e-^^ + 0(6-2-) _ ^52) 

Indeed, this result even seems to hold for complex values of A. 

I conclude with remarking that even though Takeuchi's function has been labelled a 
"Textbook Example," it provides an exciting open question for asympotic analysis. 

12 



Acknowledgements 

I thank Philippe Flajolet for bringing this problem to my attention. 



References 

[DeBruijn 1961] N. G. de Bruijn. Asymptotic Methods in Analysis, North Holland, Ams- 
terdam, 1961. 

[Knuth 1991] Donald E. Knuth. Textbook Examples of Recursion. In Artificial In- 
telligence and Theory of Computation, Academic Press, London, pages 
207-229, 1991. 

[Kuczma 1990] M. Kuczma, B. Choczewski, and R. Ger, Iterative Functional Equations, 
Encyclopedia of Mathematics and its Applications, Cambridge Univer- 
sity Press, 1990. 

[Moser 1955] Leo Moser and Max Wyman. An Asymptotic Formula for the Bell Num- 
bers. In Transactions of the Royal Society of Canada 49, 49-53, 1955. 

[Takeuchi 1978] Ikuo Takeuchi. On a Recursive Function That Does Almost Recursion 
Only. Memorandum, Musahino Electrical Communication Laboratory, 
Nippon Telephone and Telegraph Co., Tokyo, 1978. 

[Takeuchi 1979] Ikuo Takeuchi. Dai-Ni-Kai LISP Kontesuto [On the Second LISP Con- 
test]. In Joho Shori 20, pages 192-199, 1979. 



13 



